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We examine the question of whether neutralinos produced at the LHC can be 
shown to be the particles making up the astronomically observed dark matter. 
If the WIMP alllowed region lies in the SUGRA coannihilation region, then a 
strong signal for this would be the unexpected near degeneracy of the stau and 
neutralino i.e., a mass difference AM ~ (5 — 15) GeV. For the mSUGRA model 
we show such a small mass difference can be measured at the LHC using the 
signal 3r+jet+.E™ lss . Two observables, opposite sign minus like sign pairs and 
the peak of the tt mass distribution allows the simultaneous determination 
of AM to 15% and the gluino mass Mg to be 6% at the benchmark point of 
M g =850 GeV, A =0, fi >0 with 30 fb _1 . With 10 fb _1 , AM can be determined 
to 22% and one can probe the parameter space up to mi/2=700 GeV with 100 
fb- 1 . 



1. Introduction 

Supersymmetry (SUSY) offers the possibility of solving a number of the- 
oretical problems of the Standard Model (SM). Thus the cancelations im- 
plied by the bose-fermi symmetry resolves the gauge hierarchy problem, 
allowing one to consider models at energies all the way up to the GUT 
or Planck scale. Further, using the SUSY SM particle spectrum with one 
pair of Higgs doublets, (pairs of Higgs doublets are needed on theoretical 
grounds to cancel anomalies and on phenomenological grounds to give rise 
to both u and d quark masses) the renormalization group equations (RGE) 
show that grand unification of the SM gauge coupling constants occurs at 
M G ^ 10 16 GeV opening up the possibility of SUSY GUT models. How- 
ever, none of this can actually occur unless a natural way of spontaneously 
breaking SUSY occurs, and this very difficult to do with global supersym- 
metry. The problem was resolved by promoting supersymmetry to a gauge 
symmetry, supergravity (SUGRA), 1 where spontaneous breaking of super- 
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symmetry can easily occur. One can then build SUGRA GUT models 2 ' 3 
with gravity playing a key role in the construction. A positive consequence 
of this promotion was that the RGE then show that the breaking of super- 
symmetry at the GUT scale naturally leads to the required SU(2) x U(l) 
breaking at the electroweak scale, thus incorporating all the successes of 
the SM, without any prior assumptions of negative (mass) 2 terms. 

In spite of the theoretical successes of SUGRA GUTs, there has been 
no experimental evidence for its validity except for the verification of grand 
unification (which has in fact withstood the test of time for over a decade). 
However, one expects SUSY particles to be copiously produced at the LHC. 
Further, models with R parity invariance predict the lightest neutralino 
x5 to be a candidate for the astronomically observed dark matter (DM) 
and models exist 4 consistent with the amount of dark matter observed 
by WMAP, 5 and being searched for in the Milky Way by dark matter 
detectors. Thus it is possible to build models that both cover the entire 
energy range from the electroweak scale to the GUT scale and go back in 
time to ~ 1CT 7 seconds after the Big Bang when the current relic dark 
matter was created. 

The question then arises can we verify if the dark matter particles in the 
galaxy is the neutralino expected to be produced at the LHC? In principle 
this is doable. Thus assuming the DM detectors eventually detect the dark 
matter particle, they will measure the mass and cross sections, and these 
can be compared with those measured at the LHC. However, this may 
take a long time to achieve. More immediately, can we look for a signal at 
the LHC that is reasonably direct consequence of the assumption that the 
neutralino is the astronomical DM particle and in this way experimentally 
unify particle phenomena with early universe cosmology? To investigate this 
question it is necessary to choose a specific SUSY model, and for simplicity 
we consider here the minimal mSUGRA (though a similar analysis could 
be done for a wide range of other SUGRA models) . 

2. The mSUGRA model 

The mSUGRA model depends on four soft breaking parameters and one 
sign. These are; mi/ 2 (the universal gaugino soft breaking mass at Mq); 
mo (the universal scalar soft breaking mass at Mq); Aq (the universal cubic 
soft breaking mass at Mq); tan/3 =< H 2 > / < Hi > at the electroweak 
scale (where < H 2 > gives rise to u quark masses and < Hi > to d quark 
masses); and the sign of fi parameter (where fi appears in the quadratic 
part of the supcrpotential = nHiH-i). 
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Fig. 1. The feynman diagrams for annihilation of neutralino dark matter in the early 
universe 

Current experimental data significantly constrains these parameters, the 
main accelerator constraints are: The Higgs mass tojj > 114 GeV; 6 the 
lightest chargino mass M-± > 104 GeV; the b — > 57 branching ratio 2.2 
xlO -4 < Br(b — » sj) < 4.5 x 10~ 4 ; 7 and the muon g-2 anomaly 8 which 
now deviates from the SM prediction by 3.4 a. The astronomical constraint 
is the WMAP determination of the amount of dark matter and we use here 
a 2 (j range: 5 

0.094 < %o/i 2 < 0.129. (1) 

The WMAP constraint limits the parameter space to three main regions 
arising from the diagrams of Fig.l. (1) The stau- neutralino (fi — Xi) coan- 
nihilation region. Here too is small and TO1/2 < 1.5 TcV. (2)Thc focus 
region where the neutralino has a large Higgsino component. Here TO1/2 is 
small and too > 1 TcV. (3) The funnel region where annihilation proceeds 
through heavy Higgs bosons which have become relatively light. Here both 
too and TO1/2 are large. 

(In addition there is a small bulk region) Note that a key element in the 
coannihilation region is the Boltzman factor from the annihilation in the 
early universe at kT ~ 20 GeV: exp[- AM/20] where AM = M fl - M^o. 
Thus significant coannihilation occurs provided AM < 20 GeV. 
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Fig. 2. Allowed parameter space for tan/3 = 40 with Aq = and fi > 0. 




Fig. 3. SUSY production and decay channels 

The accelerator constraints further restrict the parameter space and if 
the muon g-2 anomaly maintains, fi > is prefcreed and there remains 
mainly the coannihilation region. This is illustrated in Fig. 2 which shows 
the allowed narrow coannihilation band (for the case tan/3 = 40, Ao = 0, 
fi > 0) where AM = (5 - 15) GeV and m 1/2 < 800 GeV. (There is a small 
focus region for small m 1 / 2 and mo > 1 TeV since the b — > sj constraint 
ceases to opperate at mo > 1 TeV.) 

The coannihilation band is narrow (AM = 5 — 15 GeV) due to the Boltz- 
man factor in Fig.l, the range in AM corresponding to the allowed WMAP 
range for tl^ah 2 . The dashed verticle lines are possible Higgs masses. 
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One may ask two questions; (1) Can such a small stau-neutralino mass 
difference (5-15 GeV) arise in mSUGRA, i.e. one would naturally expect 
these SUSY particles to be hundreds of GeV apart and (2) Can such a small 
mass difference be measured at the LHC? If the answers to both these ques- 
tions are affirmative, the observation of such a small mass difference would 
be a strong indication that the neutralino is the astronomical DM particle 
since it is the cosmological constraint on the amount of DM that forces the 
near mass degeneracy with the stau, and it is the accelerator constraints 
that suggests that the coannihilation region is the allowed region. 

3. Can AM be Small in SUGRA models? 

At the GUT scale mi/ 2 governs the gaugino masses, while mo the slepton 
masses. Thus, at Mq one would not expect any degeneracies between the 
two classes of particles. However, at the electroweak scale the RGE can 
modify this result. To see analytically this possibility , consider the lightest 
selectron e c which at the electroweak scale has mass 

m% =m 2 a + 0.15m 2 1/2 + (37GeV) 2 (2) 

while the \i has mass 

m| ; = 0.16m? /2 (3) 

The numerical accident that the coefficients of to^ 2 is nearly the same for 
both cases allows a near degeneracy. Thus for to = 0, the e c and Xi become 
degenerate at TO!/ 2 =(370-400) GeV. For larger TOi/ 2 , the near degeneracy 
is maintained by increasing too, so that one can get the narrow corridor in 
the too-TO!/ 2 plane seen in Fig. 2. Actually the case of the stau fi is more 
complicated since the large t-quark mass causes left-right mixing in the 
stau mass matrix and results in the fi being the lightest slepton (not the 
selectron). However, a result similar to Eqs. (1,2) occurs, with a fi — Xi 
coannihilation corridor resulting. 

We note that the results of Eqs. (1,2) depend only on the U(l) gauge 
group and so coannihilation can occur even if there were non-universal 
scalar mass soft-breaking or non-universal gaugino mass soft breaking at 
Mq. Thus, coannihilation can occur in a wide class of SUGRA models, and 
is not just a feature of mSUGRA. 
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4. Coannihilation signal at the LHC 

At the LHC, the major SUSY production processes are gluinos (g) and 
squarks (q) e.g., p + p — » g + q. These then decay into lighter SUSY parti- 
cles and Fig. 3 shows a major decay scheme. The final states involve two Xi 
giving rise to missing transverse energy E^ iss ) and four r's, two from the 
g and two from the q decay chain for the example of Fig 3. In the coanni- 
hilation region, two of the taus are high energy ( "hard" taus) coming from 
the %2 — * t ti decay (since M^o ~ 2MfJ while two are low energy ("soft" 
taus) coming from the fi — > r + \i decay since AM is small. The signal is 
thus E™ lss + jets +t's, which should be observable at the LHC detectors. 
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Number of tau pairs as a function of invariant tt mass. The difference Nos-^LS 
for mass > 100 GeV eliminating background events 



As seen above we expect two pairs of taus, each pair containing one 
soft and one hard tau from each X2 decay. Since X2 ^ s neutral, each pair 
should be of opposite sign(while SM and SUSY backgrounds, jets faking 
taus will have equal number of like sign as opposite sign events). Thus 
one can suppress backgrounds statistically by considering the number of 
oppsite sign events Nos minus the like sign events Nls< The four r final 
state has the smallest background but the acceptance and efficiency for 
reconstructing all four taus is low. Thus to implement the above ideas we 
consider here the three r final state of which two are hard and one is soft. 
(The two r final state with higher acceptance but larger backgrounds was 
discussed in, 9 and an analysis of the coannhilation signal at the ILC was 
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given in. 10 

We label three taus by their transverse energies with Ef > E% > E% 
and form the pairs 13 and 23. For signal events one of the two pairs should 
be coming from a decay and have opposite sign(OS) while the other is 
not correlated. There are two measurables that can be formed. The number 
TV and the mass of the pair M. To simulate the data we use ISAJET 7.64 11 
and PGS detector simulator. 12 Events are chosen with E™ lss and 1 jet and 
three taus with visible momenta p^ s > 40 GeV, p^ s > 40 GeV, p^ a > 20 
GeV. We assume here that it is possible to reconstruct taus with as low 
as 20 GeV. Standard Model background is reduced by requiring E^ 1 > 100 
GeV, E™ iss > 100 GeV with tevatron results, E^ 1 + E™ iss > 400 GeV. We 
also assume rate of a jet faking a r (fj^ T ) to be fj-> T = 1% (with a 20% 
error in fj^ T ) consistent with Tevatron results. 

Fig 4. shows the number of events as a function of the rr mass for 
gluino mass M s = 850 GeV and AM = 9 GeV and 20 GeV. One sees that 
the difference Nos — -Wls cancels out as expected for rr mass >100 GeV 
(consistent with the fact that the signal events are expected to lie below 
100 GeV). 

Fig. 5 shows the behavior of Nos-ls a s a function of AM and M~ g . The 
central black line is for the assumed 1%. rate for jets faking a r, the shaded 
region around it is for a 20% uncertainty in fj^ T - One sees that provided 
this uncertainty is not large, it produces only a small effect. 
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Fig. 5. Nos—LS as function of AAf (left graph) and as a function of M§ (right graph). 
The central black line assumes a 1% fake rate, the shaded area representing the 20% 
error in the fake rate 
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Figs 4 and 5 show two important features. First, Nos-ls increases 
with AM (since the t acceptance increases) and Nos-ls decreases with 
Mg (since the production cross section of gluinos and squarks decrease with 
Mg). Second, from Fig. 4 one sees that Nos-ls forms a peaked distribu- 
tion. 9,13 The ditau peak position AfP° ak increases with both AM and Mg. 
This allows us to use the two measurables Nos-ls and M^ ak to deter- 
mine both AM and M g . Fig. 6 shows this determination for the benchmark 
case of AM=9 GeV, Mg=850 GeV, A o =0 and tan/3=40. Plotted there are 
constant values of Nos-ls and constant values of M]?° ak in the AM — M g 
plane which exhibit the above dependance of these quantities on AM and 
M g . With luminosity of 30 fb _1 one determines AM and M g with the 
following accuracy: 



Fig. 7 shows how the accuracy of the measurement changes with luminosity. 
One sees that even with 10 fb _1 (which should be available at the LHC after 
about two years running)one could determine AM to within 22%, which 
should be sufficient to know whether one is in the SUGRA coannihilation 
region. 



SAM/ AM ~ 15% 



fc; SM g /M g = 6% 



(4) 




20 



AM (GeV) 



Fig. 6. Simultaneous determination of AM and M g . The three lines plot constant 
Nqs-LS an d MrT^ k (central value and la deviation) in the Mg-AM plane for the 
benchmark point of AM=9 GeV and Mg=850 GeV assuming 30 fb _1 luminosity 
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Fig. 7. Uncertainty in the determination of AM and Mg as a function of luminosity. 

5. Conclusions 

We have examined here the question of how one might show that the Xi 
particle produced at the LHC is the astronomically observed dark matter. If 
AM, the stau-ncutralino mass difference lies in the coannihilation region of 
the SUGRA mn-mi/2 plane where AM = (5—15) GeV, this would be strong 
indication that the neutralino is the dark matter particle as otherwise the 
mass difference would not naturally be so small. We saw how it was possible 
to measure such a small mass difference at the LHC for the mSUGRA model 
using a signal of S™ lss + 1 jet+3r, and simultaneously determine the gluino 
mass Mg, provided it is possible at the LHC to reconstruct taus with 
as low as 20 GeV. With 30 fb -1 one could then determine AM with 15% 
accuracy and M g with 6%, at our benchmark point of AM =9 GeV, Mg=850 
GeV, tan/3 = 40, and ^4o=0- Even with only 10fb -1 one would determine 
AM to within 25% accuracy, sufficient to learn whether the signal is in the 
coannihilation region. 

While the analysis done here was within the framework of mSUGRA, 
similar analyses can be done for other SUGRA models provided the pro- 
duction of neutralinos is not suppressed. However, the determination of M g 
does depend on the mSUGRA universality of the gaugino masses at Mq to 
relate M^o to M g . Thus a model independent method of determining M g 
would allow one to to test the question of gaugino universality. However, it 
may not be easy to directly measure M g at the LHC for high tan (3 in the 
coannihilation region due to the large number of low energy taus, and the 
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ILC would require a very high energy option to see the gluino. 

As mentioned above, one can also measure AM using the signal £?™ lss + 
2 jets+2-r. 9 This signal has higher acceptance but larger backgrounds. 
There, with 10 fb _1 one finds that one can measure AM with 18% er- 
ror at the benchmark point assuming a separate measurement of Mg with 
5% error has been made. While we have fixed our benchmark point at 
M g = 850 GeV(i.e. m 1 / 2 =360 GeV), higher gluino mass would require 
more luminosity to see the signal. One finds that with 100 fb _1 one can 
probe m 1/2 at the LHC up to ~ 700 GeV (i.e., M g up to ~ 1.6 TeV). 

Finally, it is interesting to compare with possible measurements of AM 
at the ILC. If we implement a very forward calorimeter to reduce the two 
7 background, AM can be determined with 10% error at the benchmark 
point. 10 Thus in the coannihilation region, the determination of AM at the 
LHC is not significantly worse than at the ILC. 
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